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ON 



A quantum algorithm for combinatorial search is presented that provides a simple framework for 
utilizing search heuristics. The algorithm is evaluated in a new case that is an unstructured version 
of the graph coloring problem. It performs significantly better than the direct use of quantum 
parallelism, on average, in cases corresponding to previously identified "phase transitions" in search 
difficulty. The conditions underlying this improvement are described. Much of the algorithm is 
independent of particular problem instances, making it suitable for implementation as a special 
purpose device. 
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1 Introduction 

Quantum computers |, |, @, |, [k], ||, § use quantum 
parallelism, i.e., the ability to operate simultaneously on 
a superposition of many classical states, and interference 
among different computational paths. A measurement on 
a superposition gives a definite result, with probabilities 
determined by the amplitudes of the superposition. A 
successful algorithm is one that uses these capabilities to 
arrange for a large probability of a desired result, e.g., a 
solution to a search problem. 

However, there are two major difficulties with quantum 
computers. First, they are difficult to implement [fbst . 
Second, the physical restriction to unitary linear opera- 
tions makes quantum computers difficult program effec- 
tively. An encouraging development with respect to the 
design of algorithms is a method for efficiently factoring 
integers ||| , a problem that appears to be intractable for 
classical computers. Since this relies on specific proper- 
ties of the factoring problem, the extent to which effective 
quantum algorithms can be developed for general combi- 
natorial search remains an open question. 

To help address this question empirically, it would be 
useful to have a simple framework for exploring the use 
of search heuristics with quantum algorithms. Specifi- 
cally, such a framework should allow a variety of heuristic 
search methods to be programmed while automatically 
preserving a number of desirable properties that simplify 
potential hardware implementations. This could serve to 
bridge the gap between discussions of hardware, with the 
focus usually on bit-level operations, and abstract theo- 
retical discussions, with the focus on Turing machines and 
general unitary operations. 

There are a number of desirable properties such a 
framework should provide for a general quantum search 
algorithm. First, the number of computational steps re- 
quired should be determined a priori. This avoids the 
question of when to make the final measurement, and 



could limit the difficulties of decoherence compared to 
running a program where the number of steps required 
can vary greatly from one problem instance to the next. 
Second, it is useful if much of the complexity of the algo- 
rithm can be made independent of the details of individual 
instances. This can facilitate implemention with a special 
purpose search device. Third, the framework should eas- 
ily allow applying heuristics that incorporate additional 
knowledge about the structure of particular search prob- 
lems, in a manner analogous to the heuristics used to 
dramatically improve many classical search strategies. 

Another property is given by recent studies of the na- 
ture of classes of combinatorial search problems and re- 
lates to how methods are evaluated. Most theoretical 
analyses of search algorithms focus on worst case behav- 
ior. However, in practice it is often more important to 
examine the behavior of algorithms for typical or average 
search problems. Thus even if quantum computers are not 
applicable to all combinatorial search problems, they may 
still be useful for many instances encountered in practice. 
This is an important distinction since typical instances 
of search problems are often much easier to solve than 
is suggested by worst case analyses, though even typical 
costs often grow exponentially on classical machines. 

In fact, the hard instances are not only rare but also 
concentrated near abrupt transitions in problem behavior 
analogous to physical phase transitions || [l3| [l^]. This 
result applies to many classical search methods that use 
problem structure to guide choices, but not to generate- 
and-test, where possible solutions are examined sequen- 
tially. Similarly, the most direct use of quantum paral- 
lelism for search, in which all states are simultaneously 
checked for consistency and then a measurement made, 
is equivalent to generate-and-test and does not exhibit 
the transition behavior. This limitation may apply to 
any search algorithm that uses quantum parallelism with- 
out also using interference [Q. Thus a check on whether 
a quantum algorithm is fully exploiting problem struc- 
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ture, through interference among different computational 
paths, is whether it exhibits the transition behavior. This 
is another useful property to incorporate in a general 
search framework. 

In this paper, we briefly describe a previously proposed 
general search mechanism and discuss how it meets these 
criteria as well as some implementation issues. To exam- 
ine its generality, its behavior is then evaluated empiri- 
cally for a new problem ensemble, motivated by the NP- 
complete graph coloring problem. As with other previ- 
ously studied problems, this shows a significant enhance- 
ment of probably to obtain a solution compared to the 
direct use of quantum parallelism and the transition be- 
havior. We give a more detailed look at the underlying 
basis for this improvement to see how heuristics might be 
usefully applied. Finally, some open issues are presented. 

2 A Quantum Search 

Combinatorial search can be viewed as finding, from 
among N given items, a set of size L satisfying speci- 
fied constraints. Such a set is a solution to the problem. 
The constraints, in turn, can be specified by nogoods, sets 
of items that are inconsistent. Because the supersets of 
a nogood are also nogood, sets of items can be viewed as 
forming a lattice structure consisting of levels from to 
N. Level i in this lattice contains all sets of size i, and 
these are linked to their supersets at the next level and 
subsets at the previous one. This lattice, describing the 
consistency relationships among sets, is the deep structure 
of the combinatorial search problem ^4|. This structure 
for N = 4 is shown in Fig. [l]. Solutions are found among 
the sets at level L. Notationally, we denote the size of a 
set s by \s\. 

This abstract description of search is less commonly 
used than other representations, which are more compact 
and efficient for classical search algorithms. It is intro- 
duced here as a useful basis for quantum searches and 
because it applies to many search problems. These in- 
clude constraint satisfaction problems (CSPs) [fflj such 
as graph coloring and satisfiability. For example, in col- 
oring an n-node graph with c colors an item is a pair (y, k) 
consisting of a node v in the graph and a color k for it. 
Thus there are N = nc items and a solution consists of 
L = n such items that give a unique color to each node 
and distinct colors to each pair of nodes in the graph that 
are linked by an edge. Each edge is a constraint which 
gives c nogoods, each consisting of a pair of items with 
the same color for both of the nodes linked by that edge. 
This search problem is known to be NP-complete for a 
fixed c (at least equal to 3) as n grows. From this ex- 
ample, we see that an interesting scaling regime for com- 
binatorial search is for the nogoods of the constraints to 
have a fixed, small size (e.g., at level 2 in the lattice for 




Figure 1: Structure of the set lattice for a problem with four 
items. The subsets of {1, 2, 3, 4} are grouped into levels by size 
and lines drawn between each set and its immediate supersets 
and subsets. The bottom of the lattice, level 0, represents the 
single set of size zero, the four points at level 1 represent the 
four singleton subsets, etc. 

graph coloring) while the number of items and the size of 
solutions grows linearly with problem size. 

There are many paths through the lattice from small 
sets to the larger sets at the solution level. For example, 
if N = 4 and L — 2, there are two paths from level zero to 
each set at level 2 in the lattice. E.g., the set {1, 2} is ob- 
tained via the paths 0-f{l}- 2} and 0->{2}-»{l, 2}. 
These multiple paths can be used to create interference 
with a quantum computer operating on a superposition 
of sets. One way to do this uses the fact that for search 
problems in NP it is relatively easy to test whether a 
particular set satisfies the given constraints. Thus, we 
can adjust the phase of the amplitudes along each path 
based on whether the associated set is consistent with the 
constraints. These phase adjustments attempt to create 
destructive interference among paths leading to nonsolu- 
tions and constructive interference for solutions. To the 
extent this is successful, amplitude will be concentrated 
into solution states, giving a relatively high probability to 
obtain a solution in the final measurement. A quantum 
computer with N bits can simultaneously operate on su- 
perpositions of all 2 N subsets of {1, ... , N}. Thus we can 
start with a superposition consisting of small subsets, and 
move up one level at a time to the solution level, exploring 
all possible paths 

Following this general concept, there are many possi- 
ble mappings among these sets. A particularly simple 
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method is to divide the mapping into two parts. The 
first is problem-independent and moves amplitude from 
sets of a given size to the next larger size. The second is 
an adjustment to the phases of the sets at the new level. 
A simple choice for the phase adjustment process is to 
change the sign of the amplitude of each inconsistent set 
encountered. While this is not always the optimal phase 
policy, it has the advantage of being easily computable 
since it operates independently on each set and, as shown 
below, is quite effective. 

The mapping from one level to the next is a more com- 
plex quantum operation as it involves mixing the ampli- 
tude from different states. On the other hand, this map- 
ping is independent of the details of individual problem 
instances which could simplify its implementation as a 
special purpose search device instead of a program oper- 
ating on a universal quantum computer. It is this part 
of the algorithm that gives rise to interference effects by 
mixing the contribution from different paths through the 
lattice. The operation used here is motivated by breadth- 
first classical search where each consistent set is extended 
to each possible superset at the next level. Such a map- 
ping is not unitary, though it is nearly so. Thus a simple 
choice for the quantum mapping is to use that unitary op- 
erator that is closest O] to the mapping of sets equally 
to their supersets. The resulting matrix elements have a 
simple structure, depending only on the size of the inter- 
section of the corresponding sets H] . 

2.1 The Search Method 

The search algorithm starts by evenly dividing amplitude 
among the goods at a low level of the lattice. For prob- 
lems such as graph coloring where the constraint nogoods 
involve only two items at a time, a reasonable choice is to 
start at level 2, where the number of sets is proportional 
to TV 2 . Then for each level from 2 to L — 1, we adjust 
the phases of the states depending on whether they are 
consistent and map to the next level. Let ipi be the am- 
plitude of the set s at level j after completing step j — 2 of 
the algorithm. The initial condition is ip s =1/ y N goo d s 
if the set s at level 2 is good, i.e., consistent with the 
constraints, and otherwise ips — 0. Here N goo a s is the 
number of consistent sets at level 2. 

Each step of the algorithm moves up one level giving 



assigned to the set s after testing whether it is nogood, 
and the inner sum is over all sets s of size j that have 
k items in common with r. That is, p s = 1 when s is a 
good, and otherwise p s = —I. 

After L — 2 steps we measure the state, obtaining a 
single set. This set will be a solution with probability 



(L) 



(2) 



is the 



with the sum over solution sets and p(s) 

probability to obtain the set s with the measurement of 
the final state. 



2.2 Classical Simulation 

The study of the average or typical behavior of search 
heuristics relies primarily on empirical evaluation. This 
is because the complicated conditional dependencies in 
search choices made by the heuristic often preclude a sim- 
ple theoretical analysis, although phenomenological theo- 
ries can give an approximate description of some generic 
behaviors 13, [l5| P0| . Thus as a practical matter, the 
search framework described here should allow for empir- 
ical evaluation of various heuristic methods on existing 
classical computers. 

Unfortunately, the exponential slowdown and growth in 
memory required for such a simulation severely limits the 
feasible size of the search problems. To see this, note that 
(PJ) consists of a matrix multiplication on a vector of size 



equal to the number of sets at level j of the lat- 



tice, to produce a new vector of size iVj+i. A direct eval- 
uation of this mapping requires of order NjNj+i multipli- 
cations. To reach the solution level at L = O(N) requires 
mapping at relatively high levels in the lattice where 
j = r/N with -q constant. In this case lnA^ ~ Nh(r]) 
where h{rf) = — 7711177 — (1 — 77) In (1 — 77), so the classi- 
cal computation cost scales according to e 2h ^ N , which is 
2 2N when 77 = 1/2. 

The cost of the classical simulation can be reduced sub- 
stantially (though still growing exponentially) by exploit- 
ing the map's simple structure with a recursive evaluation. 
This is done by expanding the map in (^) to include all 
2 N sets of N items and then writing (pi) as 



^ i+1) = E a * E p>*P (i) 

fc=0 \rns\=k 



where the problem-independent matrix element^] 

mapping from sets s of size j to those r of size j + 1 
that have k elements in common with s, p s is the phase 

1 These values are available on the World Wide Web Jw|, online 
appendix] . 
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where x s = Psipi^ is the original state modified by the 
choice of phases for each set, which is nonzero only for 
sets of size j. 



(3) 
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From the definition of V r jk, it follows that 

V r'ik = (M - k)V rjk + (jfe + l)y rj , fc+i 



|r'[=M-l 



(4) 



The sum in this expression is only over the | r | subsets of r 
of size | y | — 1 and so can be rapidly computed for a given 
set r. From this relation, V r j.k+i can be determined easily 
given V r jk and the values for the subsets of r with one 
element removed. Hence, from the values for V r jo we can 
readily determine the V r jk for k > iteratively starting 
from r = at the bottom of the lattice and moving up a 
level at a time until the values at level j+1 are determined. 

To evaluate V r jo define the matrix A, of size 2^ x 2^, 
with elements A rs equal to zero if \r fl s\ > 0, i.e., the 
sets have an element in common, and one otherwise. Let 
y r = A rs x s . Because x s is zero unless the set s is of 
size j, we have 

VrjO = Vr (5) 

The product Ax can be computed recursively. To see 
this consider the sets ordered by the value of the integer 
with corresponding binary representation, e.g., the sets 
without item TV come before those with N. For example, 
the sets for N = 3 are ordered as {}, {1}, {2}, {1, 2}, {3}, 
{1, 3}, {2, 3} and {1, 2, 3}. In this ordering, the matrix A 
has the recursive decomposition 



A 



A' A' 
A' 



(6) 



where A' is the same matrix but defined on subsets of 
{1, . . . , N — 1} and is the matrix all of whose entries 
equal zero. We can then compute 



Ax 



A'xW 



(7) 



where x^ 1 ) and x^ 2 ) denote, respectively, the first and 
second halves of the vector x (i.e., corresponding to sets 
without N and with N respectively). Thus the cost to 
compute Ax is C(N) = 2C{N - 1) + 0(2^) resulting in 
an overall cost of order N2 N . While still exponential, this 
improves substantially on the cost for the direct evalua- 
tion at high levels of the lattice. 

This calculation can be improved further in two ways. 
First, when computing y recursively store only the com- 
ponents for sets of size j + 1 or smaller. The values for 
larger sets are never used. Second, there is no need to 
explicitly compute and store the V r jk values individually. 
This is because we are only interested in a particular lin- 
ear combination of these values in (|l]) for sets r of size 
j + 1. Using (|]) this can be determined from a linear 
combination of values for sets of size j, and so on down to 
the bottom of the lattice. The net result of this process 
means the map of (Q) becomes a product of matrices, one 



for each step up the lattice, which in turn multiplies the 
vector Ax computed recursively as described above. 

The basic components of this computation, a recursive 
matrix combined with a simple relation to build up the 
required linear combinations, are not unique. Other pos- 
sibilities are useful to keep in mind as they may differ 
in their sensitivity to errors and also suggest a variety of 
quantum implementations. For instance, instead of the 
matrix A used above, defining B rs = (— l)' rns ' gives 



B 



B' B' 
B' -B' 



(8) 



which is unitary when normalized by dividing each ele- 
ment by 2 N / 2 . Thus z = Bx also has a simple recur- 
sive evaluation. This can be used as the basis for com- 
puting the V r jk by noting that (||) can be replaced with 
12k=o ( — 1) Vrjk — z r . Furthermore, instead of building 
up the values from subsets with (Hh, an alternate relation 



B rr 'V r ' jk — ^ ^ VrjlTllk\r\j 



(9) 



1=0 



where the sum over r' includes all sets and nikpo 

Em=0 (-lY n nikp*;m with 



p-l\/l\/a-l 
m — x j \x J \k — x 



(10) 



counting the number of sets r' with overlap k with s and m 
with r, given that s and r have overlap I and sizes <r and 
p respectively. These relations involve somewhat more 
computational operations on a classical machine than the 
method described above. However, by avoiding use of the 
nonunitary mapping between each set and its subsets, it 
may provide a method whereby quantum computers could 
exploit the simple structure of the mapping. 

3 Quantum Search Behavior 

The behavior of this search algorithm was examined 
through a classical simulation. While these results are 
limited to small search problems, they nevertheless give 
an indication of how this algorithm can dramatically in- 
crease the probability to find solutions compared to the 
direct use of quantum parallelism. As a check on the nu- 
merical errors, the norm in the final state was 1 to within 

io- 10 . 

As an example of the average behavior we consider a 
simple class of unstructured problems corresponding to 
graph coloring with 3 colors. As described above, a graph 
with n nodes has N = 3n variable-value pairs, and solu- 
tions consist of sets of L = n = N/3 of these items. The 
constraints in a graph coloring problem involve two items 
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Figure 2: Probability in sets at the solution level for two instances of unstructured problems with N = 15. On the left, the 
problem has m = 40, with 17 solutions and P so in = 0.21, on the right, m = 75, with a single solution and P so in = 0.057. The 
solutions are drawn as large gray points. The remaining sets are small black points, 
ratio Psoin 

^ 0.51. . . . • ' 
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Figure 3: Scaling of the algorithm. The left plot shows the ratio of probability to find a solution, P ao i n , with the algorithm 
to that from random selection from among the solution level sets, for L — N/3 with prespecified solution for j3 = m/N of 2 
(solid), 3 (dashed), 4 (gray), 5 (dashed gray) and 6 (thick dashed gray). This is shown on a log-scale. The right plot shows the 
scaling of P ao in on a log- log plot. Each point is the average of 1000 problem instances (except only 100 were used for N = 27). 
The points include error bars indicating the standard error of estimates of the mean, which in most cases are smaller than the 
size of the plotted points. 



at a time. Thus to randomly generate an unstructured 
version of such a problem, we select m distinct sets, each 
with two items, to be the nogoods specified by the con- 
straints. There are (^) such sets to choose from. This 
random selection of problems ignores the detailed struc- 
ture of the constraints for graph coloring, but gives a 
wider range of possible m values for the small problems 
considered here. 

Since the quantum algorithm is incomplete, i.e., it can 
find a solution but never determine that no solution ex- 
ists, we consider only problems with a solution. Soluble 
problems are rare when there are many nogoods, so for 
simplicity we use a prespecified solution, i.e., before se- 
lecting the nogoods, a particular set of size L is selected 
to be a solution. Then, the selection of m nogoods is 
only from among those that are not subsets of this spec- 
ified solution. This guarantees the problem has at least 
one solution. Qualitatively similar results are obtained 
by full random selection and then testing, via a complete 
classical search, that the problem has a solution. 



Fig. g shows how the algorithm enhances the probabil- 
ity of finding a solution for two instances: one with few 
constraints and the other with many constraints. The 
sets are ordered according to the integer whose binary 
representation corresponds to including the items in the 
set pl] |. In this case there are ( 1 5 5 ) = 3003 sets at the 
solution level, so a random selection would give a proba- 
bility of about 0.0003 to each set, much less than given to 
solutions by this algorithm. Thus the various contribu- 
tions to nonsolutions tend to cancel out among the many 
paths through the lattice. 

3.1 Scaling 

Fig. H shows the improvement over generate-and-test for 
two cases, but how well does the algorithm do over the 
ensemble of problems, and how does this behavior scale 
with increasing problem size? An appropriate choice of 
the scaling and method of generating problem instances is 
necessary for a study of average behavior so as to include 
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a significant number of hard instances. In this respect, an 
interesting scaling regime is when the number of nogoods 
at level 2 grows linearly with the size of the problem N, so 
we define (3 = m/N. This corresponds to graph coloring 
where the number of edges is proportional to the number 
of nodes, which has a high concentration of hard search 
cases ||. 

For these problems, two views of the scaling behavior 
are shown in Fig. |^. First, the linear growth on the log 
plot indicates this algorithm improves exponentially com- 
pared to the direct use of quantum parallelism. This is 
another indication of the effectiveness of this simple algo- 
rithm at concentrating amplitude into solutions. The sec- 
ond plot in the figure shows the overall probability to find 
a solution, on a log-log plot, where linear behavior cor- 
responds to a power law. Here the problem sizes feasible 
for classical simulation are too small to see the asymptotic 
behavior clearly. Nevertheless, it appears to do quite well 
for problems with few constraints. And for the remaining 
cases, a fit to a power law is closer than an exponential. 
We conclude that this algorithm improves greatly on di- 
rect use of quantum parallelism, but it is unclear whether 
that is enough to give polynomial rather than exponential 
decrease of P S oin, on average. 

3.2 Phase Transition 

Another indication of the usefulness of this quantum 
search algorithm is its behavior as the number of con- 
straints are changed for problems of a given size. This 
is shown in Fig. ^| Significantly, the figure shows this 
search algorithm also exhibits the transition behavior 
that was described above as occurring for many classi- 
cal searches ||, i.e., hard instances are concentrated near 
a change from underconstrained to overconstrained prob- 
lems. Furthermore, the expected number of trials needed 
to find a solution is essentially constant when there are 
relatively few constraints. This appears to correspond to a 
second phase transition predicted to occur among weakly 
constrained problems [^3). Thus the algorithm is using 
interference of paths to exploit problem structure in the 
same manner as sophisticated classical search methods are 
observed to do. 

This can be understood from Fig. ^, showing how the 
average value for p(s) among solutions compares with the 
expected number of solutions iV so i n for problems with dif- 
fering numbers of constraints. The average probability 
per solution grows as constraints are added. Thus the 
simple method of inverting phases at each nogood be- 
comes increasingly effective as more nogoods are added to 
the lattice. However, this growth is initially overwhelmed 
by the even faster decrease in the number of solutions, 
leading to a smaller chance to find a solution, i.e., harder 
problems. Eventually the number of solutions stops de- 
creasing so rapidly because the constraints have already 
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Figure 4: Average number of trials to find a solution, (T) = 
1/fsoin as a function of f3 for N — 15 and 24 (black and gray 
curves, respectively). Each point is the average of 1000 prob- 
lem instances, and includes error bars indicating the standard 
error of the mean, which in most cases are smaller than the 
size of the points. 

eliminated most solutions. At this point the continued 
improvement in concentration of amplitude into the re- 
maining solutions dominates, so problems get easier. Us- 
ing such a figure to display the concentration of amplitude 
into solutions can also be useful for viewing the behavior 
of different choices for the phases introduced at inconsis- 
tent sets in the lattice, perhaps in conjunction with the 
use of heuristics. 

Similar behaviors have also been seen for other prob- 
lem ensembles |lj]. These include cases correspond- 
ing to binary CSPs with two values per variable, i.e., 
L = N/2, and the NP-complete 3SAT problem where 
the constraint nogoods are sets of size 3. For these types 
of problems, slightly different ensembles can be produced 
through changes in the way problem instances are gener- 
ated. Although all these empirical observations are lim- 
ited to small problem sizes, they do suggest the search 
method applies generally to a wide range of problem en- 
sembles. 

4 Conclusion 

The lattice structure provides a general framework for 
applying quantum computers to search problems. It has 
the advantage of an a priori specification of the required 
computational steps and provides many opportunities for 
using interference among the paths through the lattice to 
each set at the solution level. While unlikely to work as 
well as special purpose algorithms for particular search 
problems, this framework may provide a basis for high 
performance search, on average, over a wide range of prob- 
lem types. 
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Figure 5: The underlying trade-off giving the transition behavior for TV = 15 and 24 (black and gray curves, respectively). The 
left plot shows the average probability in each solution set, -Psoin/A r S oin, as a function of /3. The right plot shows the average 
number of solutions N so \ n . Both plots are on a log scale. Each point is the average of 1000 problem instances, and includes 
error bars indicating the standard error of the mean, which are smaller than the size of the points. 



One important feature of this search framework is the 
ability to incorporate additional knowledge about the par- 
ticular problem structure or other search heuristics. This 
is readily included as a modification to the choice of 
phases because any such choice is guaranteed to be a uni- 
tary operation and can operate independently on each 
state. Changes to the mapping from one level to the next 
are more complicated due to the requirement to maintain 
unitarity (as well as computational simplicity). Another 
way to incorporate heuristics is by changing the initial 
conditions. In the method reported here, initially all am- 
plitude is in small consistent sets. The hope then is that 
this concentration in consistent sets is maintained as am- 
plitude is moved up through the lattice. Instead, one 
could start with amplitude spread among all sets and try 
to arrange for it to be concentrated into consistent sets 
only by the time it reaches the solution level. In some 
cases, this is observed to increase the concentration into 
solutions. 

As a possible extension to this algorithm, it would be 
interesting to see whether the nonsolution sets with rela- 
tively high probability could be useful also, e.g., as start- 
ing points for a local repair type of search |ljj. If so, the 
benefits of this algorithm would be greater than indicated 
by its direct ability to produce solution sets. This may 
also suggest similar algorithms for the related optimiza- 
tion problems where the task is to find the best solution 
according to some metric, not just one consistent with the 
problem constraints. 

There remain a number of important questions. First 
is the issue of implementation of the map from one level 
of the lattice to the next in terms of more elementary op- 
erations that are physically realizable. This could lead to 
the construction of special purpose search devices for the 
set manipulations used in the problem-independent map- 
ping between levels of the lattice. Second, how are the 
results degraded by errors and decoherence, the major dif- 
ficulties for the construction of quantum computers p6[? 



While there are some quantum approaches to error con- 
trol ||, |22| a nd studies of decoherence in the context of 
factoring U it remains to be seen how these problems 
affect the framework presented here. Third, it would be 
useful to have a theory for asymptotic behavior for large 
N, even if only approximately in the spirit of mean-field 
theories of physics. This would give a better indication of 
the scaling behavior than the classical simulations, neces- 
sarily limited to small cases, and may also suggest better 
phase choices. Considering these questions may suggest 
simple modifications to the quantum map to improve its 
robustness and scaling. There thus remain many options 
to explore for using the deep structure of combinatorial 
search problems as the basis for general quantum search 
methods. 
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